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. The e~ and e + energy spectra from the SN1987A Supernova neutrino 

' burst interactions are calculated and compared to the observed spectra in 

Kamiokande-II and 1MB experiments. Neutrino oscillations in Supernova 
and regeneration effects in the Earth, for four combinations of neutrino 
mass hierarchy (Direct/Inverted) and the value of the mixing angle #13 
(Large/Small), are taken into account. The influence of the (anti)neutrino 
production spectra in Supernova on the observed, in Kamiokande-II and 
1MB detectors, e~/e + spectra is shown. 
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1. Introduction 



When a thermal pressure of the heat, resulting from the sequence of 
nuclear fusion reactions taking place in a star, cannot longer balance the 
gravitational attraction, the last stage of a massive star evolution begins. 
For sufficiently massive stars (mass > 10M@), the collapse of their iron core 
with a mass exceeding the Chandrasekhar limit creates a type II Supernova 
(SN), which explosion is an intense source of neutrinos. It is estimated 
that, about ~99% of the SN type II binding energy (about ~10 53 erg) 
is radiated as (anti)neutrinos of all flavors. The SN neutrino flux is so 
enormous, that even (anti)neutrino signals from SN explosions outside our 
Galaxy can be detected on the Earth. Up to now, there was only one 
case of a detection of extra-solar system neutrinos, namely neutrinos from 
the SN1987A, which had exploded in the Large Magellanic Cloud. Two 
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terrestrial detectors, Kamiokande-II and 1MB [2], were able to observe 
20 events of the SN1987A neutrino interactions in total. This observation 
confirms the main features of the mechanism of the stellar evolution. 



2. The Supernova 1987A Signal 

The number of observed neutrino events, in the terrestrial detectors, 
and their energy spectrum depend on: (1) the assumed neutrino production 
fluxes, (2) the changes in neutrino fluxes due to the neutrino flavor mixing 
(oscillations) in both, the Supernova and in the Earth, and (3) the detector 
properties - including detection efficiency and the dead time. 

All these effects were taken into account in our calculations of the elec- 
tron and positron energy spectra from the SN1987A neutrino interactions 
in Kamiokande-II and 1MB detectors (these detectors were not able to dis- 
tinguish between electrons and positrons), presented in this paper. 

For the purpose of our simulations, we assume that the distance between 
the SN1987A and the Earth is 52kpc and that the total energy emitted in 
the core collapse of the star in form of (anti)neutrinos is Ltot — 3 x 10 53 erg , s. 



2.1. Primary Spectra 

We focus on the spectral characteristics of the time-integrated (anti)neu- 
trinos fluxes. The detailed spectral shape of the production fluences is not 
well known. Such fluxes are often expressed in the form of Fermi-Dirac 
distributions with the chemical potential equal to zero (see also 0), 

120 L a E 2 

e E / Ta + l' 

where a = v e , v^, v T , z/ e , u^, v T and E represents the energy of (anti)neutri- 
nos. The L a is the total energy released in various flavors of (anti)neutrinos 
(L a ~ L to4 /6, where an equipartition of the energy emitted in different 
species is assumed), the T a is the temperature of the v a gas in the (an- 
ti)neutrino sphere. In our calculations, we assume the following, "standard" , 
hierarchy of temperatures: T Ue = 3.5MeV, T Pe = 5MeV, T Ux ^ x = 8MeV, 
where u x and v x mean v^, v T and u^, v T , respectively. These temperatures 
corresponds to the following, approximate average values of the neutrino 
energies: E Ue = llMeV, E Pe = WMeV, E UxtPx = 25MeV. 

Using a nominal Fermi-Dirac distribution to approximate the spectrum 
is physically motivated because a truly thermal neutrino flux would follow 
this behavior. Numerical simulations of the Supernova explosion suggest, 
however, that the production fluences can be "pinched", meaning that their 
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Lv x j> x 


Lu x ,v x 


Accretion 


11.0 


13.5 


15.0 


4.0 


4.5 


3.0 


1.4 


1.4 


Cooling 


13.5 


16.5 


18.0 


2.7 


2.7 


2.7 


0.7 


0.7 



Table 1. The values of the parameters, of the "Analytic Fit Functions", used in 
this paper. For details see text and [S]. 



low- and high-energy parts are suppressed relative to the Fermi-Dirac spec- 
trum of the same average energy. Therefore, we consider also an alternative 
"Analytic Fit Function" 4 for which analytic simplicity is the main moti- 
vation, 

F°(E) = — • ^— + 1 ) /3a+1 . ( Ji_\ '" . , -i >'" + J t.< 



El r(/3 a + l) \E, 

where E is the energy of (anti)neutrinos, L a is the total energy released, 
E a is the average energy and @ a is the amount of spectral pinching. 

In our analysis, we sum up the spectra coming from the two phases of 
the Supernova explosion, the "accretion" phase in which (about) 10% of 
the total L tot is released, and the "cooling" phase in which the remaining 
90% of the L to t is released. The values of parameters, that we use, were 
extracted from the Garching group paper [2] and are given in Tabled 



2.2. Modifications by Oscillations 

The details of our approach can be found in so here we only recall 
the main features: (1) neutrino oscillations in SN are expressed in terms of 
flip probabilities describing their transitions in two resonance layers [Bj, (2) 
due to the large distance from the Supernova, they reach the Earth surface 
in form of incoherent mass eigenstates, (3) for the neutrino regeneration in 
the Earth, the realistic PREM I 8 Earth density profile is used. 

We consider two neutrino mass hierarchies - the Direct (Normal, mi < 
m.2 << ^3) and Inverted (7773 << mi < m<i) ones, and two values of 
013 (called "Large" and "Small", see also 0)- The following values of the 
neutrino mixing parameters were used (0, the Dirac's phase 5 = 0): 



Amli 


= 7.92 x 10~ 5 eV 2 
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= 2.4 x 10~ 3 eV 2 
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= 0.314 
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= 0.9 x 10~ 2 
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Fermi-Dirac 


Analytic Fit Function 


u e + p — > n + e + 


76%-80% 


87%-89.5% 


z/ e + -> F + e" 


6%-10% 


2%-3.5% 


i> e + O -> N + e+ 


2.5%-10% 


1.5%-2% 



Table 2. The main CC (anti)ncutrino interactions with water, together with their 
calculated percentage contributions to the total number of interactions (given 
ranges correspond to different combinations of neutrino mass hierarchy and 13 
value). 

2.3. Kamiokande-II and 1MB detectors 

Both these detectors are water Cerenkov detectors. For the description 
of the Kamiokande-II detector we use (2.14ktons of water, the efficiency 
extracted from the FIG. 3 in there, dead time negligible), while for the de- 
scription of the 1MB detector we use (6.8ktons of water, the efficiency 
extracted from the FIG.l in there, 13% dead time correction). 

The Charge Current (CC) (anti)neutrino interactions with water, which 
are considered in our calculations (together with their calculated percentage 
contributions to the total number of interactions in these detectors), are 
given in Table [21 (the cross sections for all processes come from [S] ) . 

The Fig^ shows the sketch of the positions of Kamiokande-II and 1MB 
detectors at the time of the arrival of SN1987A neutrinos. We assume 




Fig. 1. The sketch of the positions of Kamiokande-II and 1MB detectors at the time 
of the arrival of SN1987A neutrinos. 



the following distances traveled by (anti)neutrinos in the Earth during the 
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Fig. 2. The predicted e + and e~ energy spectra, in case the SN1987A production 
spectra are described by the thermal Fermi-Dirac functions, for different combi- 
nations of mass hierarchy and O13 (DL - Direct mass hierarchy and Large Q13, 
IL - Inverted mass hierarchy and Large Q13, DS - Direct mass hierarchy and 
Small 0i3, IS - Inverted mass hierarchy and Small 613). The shaded areas show 
the histograms of observed SN1987A events. For details see the description in text. 



SN1987A neutrino burst (see also JO]): 4363km for Kamiokande-II and 
8535km for 1MB. 



3. Results and discussion 

The predicted e + and e~ energy spectra, from the SN1987A neutrino 
interactions, in Kamiokande-II and 1MB detectors are presented in Figs. [2] 
and |21 respectively, for the Fermi-Dirac and the " Analytic Fit Functions" 
describing the Supernova neutrino production spectra. 

The lines correspond to the four combinations of neutrino mass hierarchy 
and 0i3 value, namely DL - Direct mass hierarchy and Large O13, IL - 
Inverted mass hierarchy and Large DS - Direct mass hierarchy and 

Small 0i3, IS - Inverted mass hierarchy and Small 613. The shaded areas 
show the histograms of the observed SN1987A neutrino events, respectively, 
12 and 8, in the Kamiokande-II and 1MB detectors. 

The predicted total numbers of registered events, in both detectors, 
are given in Table 01 The numbers in parenthesis give the predicted total 
numbers of neutrino interactions N in these detectors (all possible reactions 
taken into account, no detector efficiency and no dead time considered, see 
also 

It can be seen that none of the simulated spectra, obtained with any 
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Fig. 3. The predicted e + and e~ energy spectra, in case the SN1987A production 
spectra are described by the " Analytic Fit Functions" , for different combinations 
of mass hierarchy and 0i3. For details see the description in Fig|2|and in text. 





DL 


DS 


IL 


IS 




Fermi-Dirac 


Kamiokande-II 


22.2 


21.0 


29.3 


21.0 




(30.1) 


(29.4) 


(36.5) 


(29.4) 


1MB 


22.4 


20.3 


40.2 


20.3 




(95.6) 


(93.3) 


(116.1) 


(93.3) 




Analytic Fit Function 


Kamiokande-II 


16.0 


15.8 


20.7 


15.8 




(21.0) 


(20.8) 


(25.9) 


(20.8) 


1MB 


10.3 


10.1 


15.0 


10.1 




(66.8) 


(66.1) 


(82.5) 


(66.1) 



Table 3. The predicted total numbers of registered neutrino events in Kamiokande- 
II and 1MB detectors. The numbers in parenthesis give the predicted total numbers 
of neutrino interactions N in these detectors. For details see text. 



of the considered neutrino production energy spectra, and any of the mass 
hierarchy and G13 combinations, agrees satisfactorily with the measured dis- 
tributions. The agreement is slightly better for the 1MB spectra, especially 
for the "Analytic Fit Function" neutrino production spectra. 

Moreover, with the current neutrino mixing parameters, the different 
Earth matter effects (i.e. the different distances traveled by neutrinos in 
the Earth) cannot also explain the qualitative differences in the observed 
Kamiokande-II and 1MB spectra. In fact, as we have already noticed in [3], 
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the Earth matter effects have almost no influence on the expected total 
number of neutrino interactions N in detectors, nor on the predicted event 
spectra, except on distances below about 2500km, where some variations 
within 3.5% (1.5%), for the Fermi-Dirac ("Analytic Fit Function") neutrino 
production spectrum exist (one gets the highest numbers of events when 
neutrinos do not cross the Earth at all). 

The shape of all four (DL, IL, DS and IS) simulated e + and e~ energy 
spectra is very similar. However, the number of the expected events is 
considerably higher (and the spectrum is harder) for the IL case than for the 
remaining three, for which the distributions are practically identical (note 
also that, all these numbers scale linearly with the assumed total energy 
Ltot released in form of (anti)neutrinos, and are inversely proportional to 
the square of the assumed distance to the Supernova). 

The next generations of neutrino experiments, with much bigger detec- 
tors, and improved detection techniques, should be able to observe future 
Supernova neutrino bursts, and not only confirm the main features of the 
SN mechanism, but also to allow for more precise tests of SN models. 



REFERENCES 

[1] K.S.Hirata et al., Phys. Rev. D 38 (1988) 448-458. 
[2] C.B.Bratton et al., Phys. Rev. D 37 (1988) 3361-3363. 

[3] B.Bekman, J.H oleczek, J.Kisiel, Act. Phys. Pol. B 35 (2004) 1215 
hep-ph/0403117 . 

[4] M.T.Keil, G.G.Ra ffelt, H.T.Janka, Astrophys. J. 590 (2003) 971-991 
|astro-p h/0208035|. 

[5] G.G.Raffelt, M.T.Keil, R.Buras, H.T.Janka, M. Rampp, Contribution to Proc. 

NOON 2003 |astro-p h/0303226 . 
[6] A.S.Dighe, A.Yu.Smirnov, Phys. Rev. D 62 (2000) 033007. 
[7] G.L.Fogli, E.Lisi, A.Marrone, A.Palazzo, |hep-ph/ 0506083 
[8] A.M.Dziewonski, D.L.Andcrson, Phys. Earth Planet. Inter. 25 (1981) 297. 
[9] J.Heise (SNO collaboration), PhD Thesis, 



http://www.sno.phy.queensu.ca/sno/papers/Heise_thesis_2s_colour.pdf 



[10] C.Lunardini, A.Yu.Smirnov, hep-ph/0009356 



